Abstract. Recent measurements of quiet Sun heating 



events by Krucke r & Benz ( |1998D give strong support 
to Parker's ( 198^ ) hypothesis that small scale dissipative 



events make the main contribution to the quiet heating. 
Moreover, combining their observations with the analy- 
sis by Priest et al. ( ^000 ), it can be concluded that the 



sources driving these dissipative events are also small scale 
sources, typically of the order of (or smaller than) 2000 
km and the resolution of modern instruments. Thus arises 
the question of how these small scale events participate 
into the larger scale observable phenomena, and how the 
information about small scales can be extracted from ob- 
servations. This problem is treated in the framework of a 
' simple phenomenological model introduced in Krasnosel- 
.' skikh et al. (2001), which allows to switch between various 
^ ' small scale sources and dissipative processes. The large 
scale structure of the magnetic field is studied by means 
of Singular Value Decomposition (SVD) and a derived en- 
tropy, techniques which are readily applicable to experi- 
mental data. 
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1. Why small-scale sources ? 

The anomalously high temperature of the solar corona is 
still a puzzling problem of solar physics, despite the con- 
siderable theoretical and experimental efforts involved for 



a long time (e.g. Priest et al., |2000[ Einaudi & VeUi, |l99j. 
Since the energy release in the largest heating events 
(flares and microflares) does not supply enough power to 
heat the Corona, the statistical behavior of smaller-scale 
and less energetic but much more frequent events is an 
essential feature of the problem, as was conjectured some 
time ago by Parker (|1988D. 



An important result that supports Parker's hypoth- 
esis was reported by Krucker & Benz ( |1998| ) who have 
found from the Yohkoh / SXT observations, assuming 
that the flaring region has a constant height, that the en- 
ergy probability density has the form of a power law in 
the energy range 10^'*-10^^ ergs with an exponent about 
—2.59. Such an exponent less than —2 indicates that heat- 
ing takes place in small scales, while on the contrary an 
exponent greater than —2 indicate that large scale phe- 
nomena play a dominant role. The conclusion of Krucker 
& Benz was confirmed by Parnell & Jupp ( 200C| ), who es- 
timated the exponent to be between —2 and —2.1 making 
use of the data of TRACE and supposing that the height 
varies proportionally to the square root of the area. The 
relatively steep slope of the energy distribution found by 
these authors also suggests that the smallest flares con- 



tribute essentially to the heating. Mitra & Benz (2000) 



have discussed the same observations as Krucker & Benz 
but supposing the height variations similar to Parnell & 
Jupp, they have shown that the exponent becomes a lit- 
tle larger then previous estimate but still is smaller than 
minus two. 

Making use of the multi-wavelength analysis Benz & 



Krucker (1999) have shown that energy release mecha- 



nisms are similar in large scale loops and in the faintest 
observable events. They have also noticed that the heating 
events occur not only on the boundaries of the magnetic 
network but in the interiors of the cells also. Priest et al. 



(1998), comparing model predictions for the plasma heat- 
ing in the magnetic loop due to the distributed energy 
source with observations led to the conclusion that the 
heating is quasi-homogeneous along the magnetic loop. 
This means that the heating process does not occur in the 
close vicinity of the foot points but rather in the whole 
arc volume. Putting together these facts, it follows that 
the characteristic spatial scale of the magnetic field loops 
which supply the magnetic field dissipated is of the same 
order as the characteristic scale of the dissipation. Thus 
we may conclude that not only the dissipative process, but 
also the energy sources have small characteristic length. It 
also results from these observations that the sources are 
distributed quite homogeneously in space. 

Hence it is important to discuss the role and properties 
of sources and dissipative processes in the framework of 
simplified models. Such an approach allows to study the 
correspondence between large scale magnetic field prop- 
erties and characteristics of the small scale (eventually 
smaller than the experimental resolution) characteristics 
of sources and dissipative events, which may become a ba- 
sis of the analysis of experimental data, helping to explain 
the nature of the physical phenomena underlying the ob- 
servations. 

A phenomenological model allowing for different 
sources and physical dissipation mechanisms was proposed 
in Krasnoselskikh et al. (2001), and their effect on the 
temporal statistics of the total dissipated energy were 
studied. This model is briefly discussed in the next sec- 
tion. To study spatial properties of the magnetic fields 
and dissipative events, detailed statistical tests applicable 
both to simulations and experimental data are required. 
Tools such as the magnetic field entropy and extraction of 
the most energetic and large scale spatial/temporal eigen- 
modes by means of Singular Value Decomposition (SVD) 
are described in section III. Their application to our model 
and their ability to discriminate between various sources 
and dissipative mechanisms are discussed in section IV, 
and the final section proposes a review and a critical dis- 
cussion of the results. 
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2. Small-scale driving and dissipation 

Various phenomenological models of flare-like events and 
cooperative phenomena in the corona have been consid- 
ered i n the literature (e.g. L u fc H amilton, f 991 ; Vlahos 
et al., 1995 ; Georgoulis et al., ^001 ), mostly r elying on the 



V X (u X B) + dissipative term. 



(1) 



notion of Self-Organized Criticality (Jensen, 1998 ). Such 
models usually exhibit infinite-range spatial correlations, 
and due to the tenuosity and localization of the driving 
do not provide an appropriate framework for our purpose. 

Instead, we shall use the model introduced in (Kras- 
noselskikh et al., 2001; Podladchikova et al., 1999| ) which 
allows for a driving more distributed in space and dissi- 
pative processes relevant to heating studies. As usual, the 
model represents a simplification of the magnetohydrody- 
namic induction equation 
dB 
'dt 

The turbulent photospheric convection randomizes in 
some sense the first term of the right hand side, which 
can be replaced by various source terms with specified 
statistical properties. The dissipative terms may take into 
account different effects such as normal and anomalous 
resistivity or magnetic reconnection, which in general de- 
pend on the current density and magnetic field configu- 
ration (their meaning and differences between the two in 
this context are discussed at length in Krasnoselskikh et 
al. ( pOOl| ). 

The model is two-dimensional, the magnetic field being 
perpendicular to the grid, with periodic boundary condi- 
tions. A discrete description of the magnetic field in term 
of cells is proposed, while the currents are computed from 

j.\ B{x,y)-B{x,y + S) 

jyj 5\B{x + 5,y)~B{x,y), 

where S is cell length {6 — 1 m the following). The currents 
can be considered as propagating on the border between 
the cells, and satisfy Kirchoff 's law at each node. 

As discussed in the introduction, one may suppose that 
the source term representing the magnetic energy injection 
has a characteristic spatial comparable with that of the 
dissipative processes. Hence source terms, mimicking the 
magnetic energy injection from the turbulent photosphere, 
are assumed to have a vanishing average, and act in each 
cell of the grid at each time step. Three different types 
with different statistical properties are considered: 

— Random sources. The simplest source is given by in- 
dependent random variables 6B in the set { — 1,0,1}, 
acting in each cell. Such a source can be made dipo- 
lar by dividing the grid into two parts where random 
numbers are chosen from the similar set that have how- 
ever, positive and negative mean values, respectively, 
for each of these parts. 

— A chaotic source. Turbulence is certainly not a com- 
pletely random process, and some of its aspects are en- 
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Fig. 1. Graphical representation of the Geisel map (solid 
line). The fixed points of the map correspond to the in- 
tersections of the graph with the straight line -Bn+i = S„ 
(dashed line). 

lightened by deterministic models. The source in each 
cell evolves independently of others according to 

SB„+i = 1 - 2{SBr,)^ 

where SB G [—1; 1], which is an instance of the logistic 
(Ulam) map, well known for its chaotic properties. 
— Geisel map source. The source term may depend on 
the value of B itself. When the dissipation is absent, 
the magnetic field in each cell evolves according to a 
map 

Bn+l = f{Bn)- 

We us ed the map introduced by Geisel and Thomae, 
( 1984 ), hereafter called Geisel map (see Fig. |]), whose 
marginally stable fixed points are responsible for the 
anomalous diffusion exhibited by this map 

(B2)ocr, a<l. 

It is generally expected that magnetic field lines in 
a turbulent plasma exhibit a subdiffusive behavior, 
which is, however, more complex than described above. 

It is worth pointing out that all the source evolve in 
time in each cell independently of other cells, and that 
interaction between the fields in neighboring cells are only 
related with the dissipation effects. 

The dissipation provides the conversion of magnetic 
energy into particle acceleration and thermal energy, and 
in our model provides the coupling between the magnetic 
field elements. Dissipative processes are most important 
where a current sheet carrying strong current density has 
formed. Neglecting resistivity, which is small in the corona, 
one is left with various instabilities of magnetic field con- 
figurations which can provide the dissipation. We consider 
two of them: 
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Anomalous resistivity, which arises from the devel- 
opment of certain instabihties such as modified Bune- 
man when the electric current exceeds a certain thresh- 
old in coUisionless plasma. In our model the currents 
are simply annihilated whenever they exceed a certain 
threshold, 

Ij'I ^ Jmax- 

Reconnect ion, for which we impose which in our 
model the additional requirement that the magnetic 
field has a configuration favorable to a X-point. Hence, 
the following two conditions should be satisfied simul- 
taneously: 

\j\ = \B-B'\ > J„ax, 

B-B' < 0, 



(2) 



i -2 > i -2 



The new condition results in the existence of currents 
that can strongly overcome the critical value. 

The difference underlined here is mainly that reconnec- 
tion represents a change in equilibrium, from one topology 
(here a X-point) to another, while anomalous resistivity 
does not require any particular topology and thus may 
also act in cell interiors and not only at boundaries. An- 
other important difference is that anomalous resistivity 
provides Joule-like heating, while reconnection yields ac- 
celerated outgoing flows and thus may be associated to 
non-thermal radiation. 

When the current is annihilated, magnetic field val- 
ues in both neighboring cells, B and B' , are replaced by 
1/2{B + B'), so the density of magnetic energy dissipated 
in a single event is given by 

AE^-(B- B'f 
4 

The procedure modeling the dissipation of currents is the 
same for both anomalous resistivity and reconnection. On 
each time step, the currents satisfying the dissipation cri- 
terion are dissipated till all the currents become subcrit- 
ical (or have the same sign for the case of reconnection). 
Then, we proceed to the next time step and switch on the 
source. Indeed, dissipative processes are supposed to be 
faster than the driving ones. The total dissipated energy 
is calculated as a sum over all the dissipated currents for 
the time step considered. 

In Podladchikova et al. ( |1999D and Krasnoselskikh et 
al. (2001), the influence of the dissipative processes and 
source terms on the statistical properties of the dissipated 
energy was studied. The dissipation was shown to have 
a significant influence on the statistics of dissipated en- 
ergy. Indeed the reconnection mechanism was shown to 
yield the strongest deviation from Gaussian distribution 
in the large energies. However, the probability density of 
the dissipated energy was shown to be rather insensitive to 
the nature of the magnetic field sources. In this paper we 
would like to explore further the dependence of the large 
scale magnetic field statistical properties upon the physi- 
cal characteristics of the source and dissipation processes 



in the framework of our model. The aim of this work is 
to investigate the possibility of getting information about 
small scale magnetic energy sources making use the large 
scale magnetic field. In order to do so, we study hereafter 
the spatial complexity of the large scale field by means 
of spatial correlations, entropy and eigenmodes (Singular 
Value Decomposition). 

3. Characterization of spatial complexity 

Spatial complexity can be characterized in many different 



ways (e.g. Grassberger, 1986). Linear properties are tradi- 
tionally studied by considering the time averaged spatial 
correlation function 

C(r) = (B(x, t) B(x + r, i))x,t/(B(x, tf)^,t , (3) 
where the average is carried out over different positions 
and times (or events). We have computed the character- 
istic decay length of this correlation function for various 
sources, dissipation mechanisms, and thresholds. 

A different approach, which is used in the framework of 
image processing, is based on the Singular Value Decom- 
position (SV D) or Karhuncn-Loeve Transform, see Golub 
& van Loan ( 1996 ). For each time step, the bivariate mag- 
netic field intensity B{x,y) can be viewed as 2D image. 
We decompose this image into a set of separable spatial 
modes 

N 

B{x,y)=^Hk fk{x) gliy) . (4) 

k=l 

By making these modes orthogonal {fkf*) ~ {9kg*) = 
Sk^i, the decomposition becomes unique. The weights fik of 
these modes, also called singular values, are conventionally 
sorted in decreasing order, and are invariant with respect 
to all orthogonal transformations of the matrix B(x, y). In 
our case, the number TV of modes is equal to the spatial 
grid size. 

A key property of the SVD is that it captures large- 
scale structures in heavily weighted modes, whereas pat- 
terns that are little correlated in space are deferred to 
modes with small weights. The distribution of the singular 
values is therefore indicative of the spatial disorder: a flat 
distribution means that there is no characteristic spatial 
scale and hence, the magnetic field should not show large- 
scale patterns. Conversely, a peaked distribution s uggest s 



that there are coherent structures (Dudok de Wit, 1995 ) 
It must be stressed that this approach is, like the previous 
one, based on second order moments only, since the spatial 
modes and their singular values issue from the eigenstruc- 
ture of the spatial correlation matrix of the magnetic field. 

From the SVD modes of the 2D magnetic field, one 
can define a measure of spatial complexity, which is based 



on the SVD entropy (Aubry |1991| ). Let Eu = ^^l/Y.^l4 
be the fractional amount of energy which is contained in 
the fc'th mode. The SVD entropy can then be defined as 
1 ^ 

Y^EklogEk. (5) 



H 



lim 



logA^ 



k=l 
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Fig. 2. Correlation length dependence on the threshold of 
dissipation for both dissipation rules. Anomalous resistiv- 
ity is marked by triangles, reconnection by squares. 

The maximum value H ~ 1 is reached when spatial disor- 
der is maximum, that is when Ek = 1/N for all k. H ^ 
means that all the variance is contained in a single mode. 

The SVD can also be used as a linear filter to extract 
large scale patterns from a background with small scale 
fluctuations. To do so, one should perform the SVD and 
then in eq. ^ sum over the strongest modes only, to obtain 
a filtered magnetic field. There is obviously some arbitrari- 
ness involved in the identification of what we call strong 
modes, but the process can be automated by usin g robu st 
selection criteria, see for example Dudok de Wit (1995). 




Fig. 3. Averaged magnetic field spatial correlation coeffi- 
cient in log-linear plot (solid lines). Dotted lines represent 
the best fits by exponential function. The results are ob- 
tained for random and subdiffusive sources with a thresh- 
old of dissipation jmax = 1. a) random source, anomalous 
resistivity dissipation, correlation length L =^ 17; b) ran- 
dom source, reconnection dissipation, correlation length 
L — 19; c) subdiffusive source, anomalous resistivity dis- 
sipation, correlation length L — 23; d) subdiffusive source, 
reconnection dissipation. 



4. Spatial complexity and properties of the source 
and dissipation 

4-.1. Spatial correlations 

Spatial correlation function as defined by Eq. ^ were av- 
eraged in time after the system has reached a stationary 
state. 

For the small grids, of the order of 30 x 30, the cor- 
relation function decays as a power-law. It was shown in 



Krasnoselskikh et al. 2001 the probability density of the 



total dissipated energy also decays as a power-law in this 
case. However the resemblance to a self organized critical 
behavior is only an artifact of the small grid size, when the 
correlation length is comparable with the grid size, and 
disappears for large enough grids. Indeed, for grid sizes 
around 100 x 100 (or greater), the time averaged corre- 
lations functions have an almost exponential decay, while 
the dissipated energy has a quasi-Gaussian distribution. 

In case of exponential correlation functions, a correla- 
tion length is easily defined as the parameter L such that 
C(r) = exp(-r/i). 

We have found in our simulations that the correlation 
length L remains almost constant when the grid size is in- 



creased above 200 x 200, and furthermore that L is much 
smaller than the linear grid size. In such a case, it is legit- 
imate to expect that results do not depend significantly 
on the grid size or boundary conditions. In the remainder 
of this paper all results are presented for a grid 400 x 400 
and a small threshold is used {jmax = 1, which is of the 
order of yJ{5B'^)). 

As shown on Fig. |^, for the random magnetic field 
source the average correlation length is larger for the re- 
connection type dissipation than for anomalous resistiv- 
ity. This can be explained by the presence of supercritical 
currents > jmax which may exist because of condition (H) . 
Moreover in both cases the correlation lengths are decreas- 
ing functions of jmax (see Fig. ^) . These results are hardly 
changed if the random source is replaced by the Ulam 
map source (see previous section). Only in the case of a 
magnetic field following the Geisel map, and when dissipa- 
tive processes are ruled by anomalous resistivity (Fig. ||c), 
the average correlation length is a little bit larger than 
in the previous cases (see Fig. ^c). When dissipation oc- 
curs through reconnection, the average correlation func- 
tion seems featureless (see Fig. |3|i) (to the least neither 
exponential nor power-law), and no correlation length can 
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Fig. 4. Typical magnetic field for random source and reconnection, with jmax = 1. a) Magnetic field a.t t — 1000, with 
entropy H = 0.79; b) Zoom of the precedent image; c,d) filtered magnetic field, with 20 eigenmodes, a,t t — 1000 and 
t = 2000 {H = 0.807). 



be unambiguously computed for comparison with previous 
results. 

Thus the correlation function seems to have difficul- 
ties here in indicating significant differences between the 
different processes. However, strong differences are seen 
simply by visible inspection of the magnetic field between 
for example the random and Geisel map sources (Figs. 



^ and Figs. ^). Hereafter we shall present an alterna- 
tive analysis using the spatial entropy defined in terms of 
SVD eigenvalues. 
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Fig. 5. Magnetic field for the subdiffusive source and reconnection, with jmax = 1- a) Magnetic field in a transient 
state, H = 0.73; b), Magnetic field in a stationary state, H = 0.51 c) 20 modes of the precedent image; d), magnetic 
field at t = 3000, H = 0.527 



4-2. Singular values and coherent spatial modes 

As discussed in the previous section, the Singular Value 
Decomposition provides an orthogonal decomposition 
which allows to extract coherent patterns eventually ex- 
isting in the bivariate magnetic field at a given time. 

In our case, the distribution of singular values ap- 
pears to be significantly peaked (Fig. ^), showing that 
the bi-dimensional wavefield is dominated by a few spa- 
tial modes. For instance, the most energetic mode {fi{x), 
in the notation of Eq. ^ obtained from SVD of the mag- 
netic field for a Geisel source and dissipation by reconnec- 
tion clearly corresponds to a large scale coherent structure 
(Fig.0). 



The fact that the most heavily weighted modes cor- 
respond to large scale magnetic field structures is further 
illustrated by Figs. ^ andJM by comparison of the mag- 
netic field (Figs. ^3 and gp) with an approximation of 
the magnetic field containing only 20 modes (Figs. ^ and 
which has a very similar large scale structure. In both 
cases, the retained modes correspond to the fast exponen- 
tial decay of the strongest singular values, down to the 
turning point where the decay becomes slower (as in Fig. 

and the following singular values are filtered out and 
set to before the inverse SVD is performed. In particular, 
it is clear that the large scale structure of magnetic field 
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Fig. 6. Energy of the spatial eigenmod.es of magnetic field 
formed by subdiffusive source and reconnection dissipa- 
tion. 
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Fig. 7. Most energetic spatial mode /i (x) obtained by the 
same Singular Value Decomposition as in Fig. ^. 



appearing with Geisel map sources, and not with random 
sources, is retained by SVD (Fig. ^). 

However, this analysis provides only a decomposition 
of the magnetic field at a given time, and no information 
about the lifetime of these structures which is quite cru- 
cial though. Monitoring the time evolution of the system, 
it appears that the heavily weighted modes appear to per- 
sist for long times, as can be seen comparing the filtered 
magnetic field at two times separated by 2000 time steps 
(Figs He and d). Actually, it is seen on Figs |^a,b,c how 
these structures grow from the initially disordered state. 

Thus the coherent structures extracted by SVD have 
a long lifetime and produce a slow decay of the temporal 
autocorrelation function, defined by 

C(t) = {B{x, y, t + t)B{x, y, t)},,^,,*, 

as shown in Fig. |^. While small-scale structures rapidly 
appear and disappear, the large scale ones evolve slowly. 
In that sense, they are truely coherent structures. 

4-.3. Magnetic field entropy 

Quantitatively, the coherence degree of the magnetic field 
can be measured by the spatial entropy defined by Eq. ^ 
from the singular values. This definition involves a limit 
A'' — > oo, but in practice, for large enough grid sizes, it 
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Fig. 8. Averaged temporal correlation function of the 
magnetic field (solid lines) in log-linear plot calculated 
from 4 X lO'' times steps. Dotted lines represent the best 
fits by an exponential function. The results are obtained 
for random and subdiffusive sources, reconnection dissipa- 
tion, with a threshold of dissipation jmax = 1. a) random 
source, correlation time r = 202; b) the same as previous, 
but only for the first 20 modes, r = 220; c) subdiffusive 
source, r — 1435; d) the same as (c), but only for first 20 
modes, r = 2958. 



can be checked that the quantity computed for a, N x N 
subset of B 

1 ^ 



log A 



fc=i 



converges toward a well-defined limit as A^ increases. Com- 
puting the entropy Hj^ for increasing A^, we have obtained 
the curves displayed on Fig. || which show that the en- 
tropy already converges for matrix sizes about 100 x 100, 
although it seems that the convergence is faster for the 
subdiffusive source than for the random source. Thus we 
may conclude that, provided large enough grids are con- 
sidered, the entropy H is fairly independent of the grid 
size. 

The entropy has a monotonous decay in time and con- 
verges toward a finite value in the steady state (see Tab. 
0), indicating the simultaneous decrease of spatial com- 
plexity and the formation of slowly evolving large scale 
magnetic field structures. 

The major result here is that the value toward which 
the entropy converges in time exhibit significant differ- 
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Fig. 9. Partial entropies as a function of the grid size 
N. The continuous line is for the subdiffusive source, and 
the dashed one for the random source. This entropy is nor- 
malized so that H — I corresponds to maximum disorder. 



t 


100 


500 


20000 


30000 


H 


0.73 


0.69 


0.51 


0.527 



Table 1. Variation of the entropy in time, for the subdif- 
fusive source and reconnection (see also Fig. ||). 



source type 


H 


random 


0.8 


Ulam 


0.78 


Geisel 


0.53 


in the steady 


state 



Table 2. Enti 
types, and dissipation by reconnection. 



ences between the different sources are used, as summa- 
rized in Table |2[ 

5. Conclusion and discussion 

To study coronal heating due to dissipation of small-scale 
current layers, we have performed a statistical analysis 
of a simple mod el. Th e model was introduced in Kras- 
noselskikh et al. 
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and its principal difference with 
previous ones is that the system is driven by small scale 
homogeneously distributed sources acting on the entire 
grid for each time step. The idea to consider small scale 
sources is motivated by observations by Benz & Krucker 
( |1998| , |1999D that the heating occurs on the level of the 
chromosphere, thus the magnetic field structures, dissipa- 
tion of which supplies the energy for the heating, are also 
of a small scale. 

The question addressed in this paper was the following: 
if the actual measurements cannot resolve the character- 
istic scale of the heating, in what sense are the "macro- 
scopic" observable properties influenced by the properties 
of the smaller scale sources ? 



To this purpose we have carried out the comparative 
analysis of statistical estimations the large scale spatial 
characteristics of the magnetic field such as the correla- 
tion length, entropy and most energetic eigenmodes for 
the different source types that were used in the model 
(random, chaotic and intermittent with anomalous tem- 
poral diffusion). 

The "noisy" small scales were filtered out in order to 
study large scale of the magnetic field. For this purpose 
we have reconstructed the magnetic field from eigenmodes 
given by SVD that corresponds to most energetic coherent 
structures. The less energetic modes that corresponds to 
noise level were truncated. 

The results can be summarized as follows: 

The large scale spatial characteristics of the magnetic 
field such as the correlation length, entropy and most en- 
ergetic eigenmodes depend significantly on both the sta- 
tistical properties of small scale magnetic field sources and 
the dissipation mechanisms. 

— It was found that the temporal average of the correla- 
tion function is exponential, i.e the correlation length 
is finite and not infinite as supposed to be in SOC 
systems. This length is a little bit larger for the recon- 
nection dissipation and it depends on the dissipation 
threshold also. 

— With the subdiffusive (Geisel) source and reconnection 
dissipation the correlation significantly departs from 
the exponential. 

— The Singular Value Decomposition (SVD) allows to 
extract the most energetic magnetic field structures, 
which are essentially larger than the source size and 
persist for long times, supporting the idea that the 
plasma can organize on large scales while being driven 
by small scale sources. 

— Moreover, the entropy computed from the singular val- 
ues of the magnetic field generated by intermittent 
sources was found to be much smaller (about 20-30%) 
for the subdiffusive source than for other sources. The 
most intensive in space and long life structures are es- 
sentially larger in this case also. That indicate a higher 
level of organization in the system than in the random 
source case. 

The clear difference of the characteristics of spatial 
complexity in the case of Geisel map sources can be ex- 
plained in the following way. This deterministic map pro- 
duces in each cell a random-like diffusion slower than usual 
(subdiffusion) of magnetic field intensity. On the other 
hand, the dissipation produces a normal diffusion of the 
field, i.e. faster magnetic field relaxation along the spatial 
grid (on average), and relates the temporal properties of 
the source to spatial properties. This explains why sources 
with slower diffusion (Geisel) tend to form larger scale and 
longer lived structures than sources with normal diffusion 
(random, Ulam). 
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Thus wc have shown in the framework of our model 
that the large scale spatial structure of the magnetic field 
in the solar atmosphere also contains important statistical 
information about the mechanisms of the coronal heating. 
Such an information can be extracted by SVD-based tech- 
niques, which are readily applicable to experimental data 
and can be used in complement to the usual analysis of 
radiated energy. 
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